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Abstract 

Inclusive transverse momentum spectra of charged particles in photoproduction events 
in the laboratory pseudorapidity range —1.2 < rj < 1.4 have been measured up to px = 
8 GeV using the ZEUS detector. Diffractive and non-diffractive reactions have been 
selected with an average "yp centre of mass (cm.) energy of (W) = 180GeV. For diffractive 
reactions, the px spectra of the photon dissociation events have been measured in two 
intervals of the dissociated photon mass with mean values (Mx) = 5 GeV and 10 GeV. 
The inclusive transverse momentum spectra fall exponentially in the low pt region. The 
non-diffractive data show a pronounced high px tail departing from the exponential shape. 
The pt distributions are compared to lower energy photoproduction data and to hadron- 
hadron collisions at a similar cm. energy. The data are also compared to the results of a 
next-to-leading order QCD calculation. 
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1 Introduction 



HERA is a prodigious source of quasi-real photons from reactions where the electron is scattered 
at very small angles. This permits the study of photoproduction reactions at photon-proton 
centre of mass (cm.) energies an order of magnitude larger than in previous fixed target 
experiments. 

The majority of the 7p collisions are due to interactions of the proton with the hadronic 
structure of the photon, a process that has been successfully described by the vector meson 
dominance model (VDM)0. Here, the photon is pictured to fluctuate into a virtual vector 
meson that subsequently collides with the proton. Such collisions exhibit the phenomenological 
characteristics of hadron-hadron interactions. In particular they can proceed via diffractive 
or non-diffr active channels. The diffractive interactions are characterized by very small four 
momentum transfers and no colour exchange between the colliding particles leading to final 
states where the colliding particles appear either intact or as more massive dissociated states. 
However, it has been previously demonstrated that photoproduction collisions at high transverse 
momentum cannot be described solely in terms of the fluctuation of the photon into a hadron- 
like state ||, [|. The deviations come from contributions of two additional processes called direct 
and anomalous. In the former process the photon couples directly to the charged partons inside 
the proton. The anomalous component corresponds to the process where the photon couples 
to a qq pair without forming a bound state. The interactions of the photon via the hadron-like 
state and the anomalous component are referred to as resolved photoproduction, since both of 
them can be described in terms of the partonic structure of the photon H. 
In this paper we present the measurement of the transverse momentum spectra of charged 
particles produced in photoproduction reactions at an average cm. energy of (W) = 180 GeV 
and in the laboratory pseudorapidity range —1.2 < r\ < 1.4 []. This range approximately 
corresponds to the cm. pseudorapidity interval of 0.8 < r] c _ m _ < 3.4, where the direction is 
defined such that positive r/ c . m . values correspond to the photon fragmentation region. The 
transverse momentum distributions of charged particles are studied for non-diffr active and 
diffractive reactions separately. The pt spectrum from non-diffr active events is compared to 
low energy photoproduction data and to hadron-hadron collisions at a similar cm. energy. In 
the region of high transverse momenta we compare the data to the predictions of a next-to- 
leading order QCD calculation. 

The diffractive reaction (p/p — > Xp), where X results from the dissociation of the photon, was 
previously measured by the E612 Fermilab experiment at much lower cm. energies, 11.8 < 
W < 16.6 GeV ||. It was demonstrated that the properties of the diffractive excitation of the 
photon resemble diffraction of hadrons in terms of the distribution of the dissociated mass, the 
distribution of the four-momentum transfer between the colliding objects || and the ratio of 
the diffractive cross section to the total cross section M. The hadronization of diffractively 
dissociated photons has not yet been systematically studied. In this analysis we present the 
measurement of inclusive pt spectra in two intervals of the dissociated photon mass with mean 
values of (M x ) = 5 GeV and (M x ) = 10 GeV. 



1 Pseudorapidity r\ is calculated from the relation r\ = —ln(tan{9 / '2)) , where 9 is a polar angle calculated 
with respect to the proton beam direction. 
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2 Experimental setup 



The analysis is based on data collected with the ZEUS detector in 1993, corresponding to an 
integrated luminosity of 0.40 ptT 1 . The HERA machine was operating at an electron energy of 
26.7 GeV and a proton energy of 820 GeV, with 84 colliding bunches. In addition 10 electron 
and 6 proton bunches were left unpaired for background studies (pilot bunches). 
A detailed description of the ZEUS detector may be found elsewhere [[7| ||. Here, only a brief 
description of the detector components used for this analysis is given. Throughout this paper 
the standard ZEUS right-handed coordinate system is used, which has its origin at the nominal 
interaction point. The positive Z-axis points in the direction of the proton beam, called the 
forward direction, and X points towards the centre of the HERA ring. 

Charged particles created in ep collisions are tracked by the inner tracking detectors which 
operate in a magnetic field of 1.43 T provided by a thin superconducting solenoid. Immediately 
surrounding the beampipe is the vertex detector (VXD), a cylindrical drift chamber which 
consists of 120 radial cells, each with 12 sense wires running parallel to the beam axis ]9[]. 
The achieved resolution is 50 fim in the central region of a cell and 150 /zm near the edges. 
Surrounding the VXD is the central tracking detector (CTD) which consists of 72 cylindrical 
drift chamber layers organized in 9 superlayers [[K]]. These superlayers alternate between those 
with wires parallel to the collision axis and those with wires inclined at a small angle to provide 
a stereo view. The magnetic field is significantly inhomogeneous towards the ends of the CTD 
thus complicating the electron drift. With the present understanding of the chamber, a spatial 
resolution of ~ 260 /im has been achieved. The hit efficiency of the chamber is greater than 
95%. 

In events with charged particle tracks, using the combined data from both chambers, the 
position resolution of the reconstructed primary vertex are 0.6 cm in the Z direction and 0.1 cm 
in the XY plane. The resolution in transverse momentum for full length tracks is <J Pt /pt < 
0.005 • pt © 0.016 (pt in GeV). The description of the track and the vertex reconstruction 
algorithms may be found in |ll|] and references therein. 

The solenoid is surrounded by the high resolution uranium-scintillator calorimeter (CAL) di- 
vided into the forward (FCAL), barrel (BCAL) and rear (RCAL) parts |L2| . Holes of 20 x 20 cm 2 



in the centre of FCAL and RCAL are required to accommodate the HERA beam pipe. Each 
of the calorimeter parts is subdivided into towers which in turn are segmented longitudinally 
into electromagnetic (EMC) and hadronic (HAC) sections. These sections are further subdi- 
vided into cells, which are read out by two photomultiplier tubes. Under test beam condi- 



tions, an energy resolution of the calorimeter of (Je/E = 0.18/ J E( GeV) for electrons and 



<te/E = 0.35/ yE( GeV) for hadrons was measured. In the analysis presented here CAL cells 
with an EMC (HAC) energy below 60 MeV (110 MeV) are excluded to minimize the effect of 
calorimeter noise. This noise is dominated by uranium activity and has an r.m.s. value below 
19 MeV for EMC cells and below 30 MeV for HAC cells. 

The luminosity detector (LUMI) measures the rate of the Bethe-Heitler process ep — * e^p. The 
detector consists of two lead-scintillator sandwich calorimeters installed in the HERA tunnel 
and is designed to detect electrons scattered at very small angles and photons emitted along the 
electron beam direction JTB|. Signals in the LUMI electron calorimeter are used to tag events 
and to measure the energy of the interacting photon, E^, from E 1 = E e — E' e = 26.7 GeV — Eg, 
where E' e is the energy measured in the LUMI. 
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3 Trigger 



The events used in the following analysis were collected using a trigger requiring a coincidence 
of the signals in the LUMI electron calorimeter and in the central calorimeter. The small 
angular acceptance of the LUMI electron calorimeter implied that in all the triggered events 
the virtuality of the exchanged photon was between 4 • 1CT 8 < Q 2 < 0.02 GeV 2 . 
The central calorimeter trigger required an energy deposit in the RCAL EMC section of more 
than 464 MeV (excluding the towers immediately adjacent to the beam pipe) or 1250 MeV 
(including those towers). In addition we also used the events triggered by an energy in the 
BCAL EMC section exceeding 3400 MeV. At the trigger level the energy was calculated using 
only towers with more than 464 MeV of deposited energy. 

4 Event selection 

In the offline analysis the energy of the scattered electron detected in the LUMI calorimeter 
was required to satisfy 15.2 < E' e < 18.2 GeV, limiting the 7p cm. energy to the interval 
167 < W < 194 GeV. The longitudinal vertex position determined from tracks was required 
to be —35 cm < Z vertex < 25 cm. The vertex cut removed a substantial part of the beam gas 
background and limited the data sample to the region of uniform detector acceptance. The 
cosmic ray background was suppressed by requiring the transverse momentum imbalance of the 
deposits in the main calorimeter, Passing, relative to the square root of the total transverse 
energy, y/E^, to be small: Passing I 'VEt < 2V GeV. 




The data sample was divided into a diffractive and a non-diffractive subset according to the 
pseudorapidity, i] max , of the most forward energy deposit in the FCAL with energy above 
400 MeV. The requirement of ri max < 2 selects events with a pronounced rapidity gap that are 
predominantly due to diffractive processes (~ 96% according to Monte Carlo (MC) simulation, 
see section 6). The events with r] max > 2 are almost exclusively (ss 95%) due to non-diffractive 
reactions. The final non-diffractive data sample consisted of 149500 events. 
For the diffractive data sample (i] max < 2) an additional cut i] max > —2 was applied to suppress 
the production of light vector mesons V in the diffractive reactions 7J9 — > Vp and 7p — > 
VN, where N denotes a nucleonic system resulting from the dissociation of the proton. The 
remaining sample was analyzed as a function of the mass of the dissociated system reconstructed 
from the empirical relationship 



The above formula exploits the fact that in tagged photoproduction the diffractively excited 
photon state has a relatively small transverse momentum. The total hadronic energy, E, 
and longitudinal momentum Pz = E ■ cosO were measured with the uranium calorimeter by 
summing over all the energy deposits of at least 160 MeV. The correction factors A = 1.7 
and B = 1.0 GeV compensate for the effects of energy loss in the inactive material, beam pipe 
holes, and calorimeter cells that failed the energy threshold cuts. The formula was optimized to 
give the best approximation of the true invariant mass in diffractive photon dissociation events 
obtained from MC simulations, while being insensitive to the calorimeter noise. The diffractive 
data were analyzed in two intervals of the reconstructed mass, namely 4 < Mx rec < 7 GeV and 
8 < M x rec < 13 GeV. According to the MC simulation the first cut selects events generated 
with a mass having a mean value and spread of (M x gen) = 5 GeV and r.m.s. = 1.8 GeV. The 




Mxrec ^A-^E 2 -P 2 + B = A-^/(E + P z )-E 1 + B. 
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second cut results in (M X gen) = 10 GeV and r.m.s. = 2.3 GeV. Details of the MC simulation 
are given in section [| The final data sample consisted of 5123 events in the lower Mx interval 
and of 2870 events in the upper interval. 

The contamination of the final data samples from e-gas background ranges from < 0.1% (non- 
diffractive sample) to « 10% (diffractive sample, (Mx) = 5 GeV). The p-gas contribution is 
between 1% (non- diffractive sample) and 2% (diffractive sample, (Mx) = 5 GeV). The e-gas 
background was statistically subtracted using the electron pilot bunches. A similar method 
was used to correct for the p-gas background that survived the selection cuts because of an 
accidental coincidence with an electron bremsstrahlung (ep — > jep). A large fraction of these 
background events were identified using the LUMI detector, since the energy deposits in the 
electron and photon calorimeters summed up to the electron beam energy. The identified 
background events were included with negative weights into all of the distributions in order to 
compensate for the unidentified part of the coincidence background. A detailed description of 
the statistical background subtraction method may be found in [EL 111. 



5 Track selection 

The charged tracks used for this analysis were selected with the following criteria: 

• only tracks accepted by an event vertex fit were selected. This eliminated most of the 
tracks that came from secondary interactions and decays of short lived particles; 

• tracks must have hits in each of the first 5 superlayers of the CTD. This requirement 
ensures that only long, well reconstructed tracks are used for the analysis; 

• — 1.2 < rj < 1.4 and px > 0.3 GeV. These two cuts select the region of high acceptance 
of the CTD where the detector response and systematics are best understood. 

Using Monte Carlo events, we estimated that the efficiency of the charged track reconstruction 
convoluted with the acceptance of the selection cuts is about 90% and is uniform in p?. The 
contamination of the final sample from secondary interaction tracks, products of decays of short 
lived particles, and from spurious tracks (artifacts of the reconstruction algorithm) ranges from 
5% at pt = 0.3 GeV to 3% for p? > 1 GeV. The inefficiency and remaining contamination of 
the final track sample is accounted for by the acceptance correction described in the following 
section. 

The transverse momenta of the measured tracks displayed no correlation with rj over the con- 
sidered interval and were symmetric with respect to the charge assigned to the track. 



6 Monte Carlo models 

For the acceptance correction and selection cut validation we used Monte Carlo events gen- 
erated with a variety of programs. Soft, non-diffractive collisions of the proton with a VDM 



type photon were generated using HERWIG 5.7 with the minimum bias option [15]. The gen- 
erator was tuned to fit the ZEUS data on charged particle multiplicity and transverse energy 
distributions. For the evaluation of the model dependence of our measurements we also used 



events from the PYTHIA generator with the soft hadronic interaction option |Tg]. Hard re- 
solved and direct subprocesses were simulated using the standard HERWIG 5.7 generator with 
the lower cut-off on the transverse momentum of the final-state partons, pTmin, chosen to be 
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2.5 GeV. For the parton densities of the colliding particles, the GRV-LO |17| (for the photon) 
and MRSD'_ |L8| (for the proton) parametrisations were used. As a cross-check we also used 



hard "fp scattering events generated by PYTHIA with pTmin — 5 GeV. The soft and hard MC 
components were combined in a ratio that gave the best description of the transverse momen- 
tum distribution of the track with the largest pr in each event. For prmin = 2.5 GeV, the hard 
component comprises 11% of the non-diffractive sample and for PTmin — 5 GeV only about 3%. 
Each diffractive subprocess was generated separately. The diffractive production of vector 
mesons (p, u, <ft) was simulated with PYTHIA. The same program was used to simulate the dou- 
ble diffractive dissociation ('yp — > XN). The diffractive excitation of the photon (-jp — > Xp) was 
generated with the EPDIF program which models the diffractive system as a quark-antiquark 
pair produced along the collision axis Final state QCD radiation and hadronization were 
simulated using JETSET [[HJ. For the study of systematic uncertainties, a similar sample of 
events was obtained by enriching the standard PYTHIA diffractive events with the hard com- 
ponent simulated using the POMPYT Monte Carlo program (hard, gluonic pomeron with the 
direct photon option) [^0] . 

The MC samples corresponding to the diffractive subprocesses were combined with the non- 
diffractive component in the proportions given by the ZEUS measurement of the partial pho- 
toproduction cross sections |J. The MC events were generated without electroweak radiative 
corrections. In the considered W range, the QED radiation effects result in « 2% change in the 
number of measured events so that the effect on the results of this analysis are negligible. The 
generated events were processed through the detector and trigger simulation programs and run 
through the standard ZEUS reconstruction chain. 



7 Acceptance correction 

The acceptance corrected transverse momentum spectrum was derived from the reconstructed 
spectrum of charged tracks, by means of a multiplicative correction factor, calculated using 
Monte Carlo techniques: 

C(pr) = (t^— ■ ^ L )/(lJ— ■ ^% 

^ gen ev U-PTgen ^recev (^PTrec 

N gen denotes the number of primary charged particles generated with a transverse momentum 
Prgen in the considered pseudorapidity interval and N genev is the number of generated events. 
Only the events corresponding to the appropriate type of interaction were included, e.g. for the 
lower invariant mass interval of the diffractive sample only the Monte Carlo events correspond- 
ing to diffractive photon dissociation with the generated invariant mass 4 < M Xge n < 7 GeV 
were used. N rec is the number of reconstructed tracks passing the experimental cuts with a 
reconstructed transverse momentum of prrec, while N recev denotes the number of events used. 
Only the events passing the trigger simulation and the experimental event selection criteria 
were included in the calculation. To account for the contribution of all the subprocesses, the 
combination of the MC samples described in section |6] was used. 
This method corrects for the following effects in the data: 

• the limited trigger acceptance; 

• the inefficiencies of the event selection cuts, in particular the contamination of the diffrac- 
tive spectra from non-diffractive processes and the events with a dissociated mass that 
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was incorrectly reconstructed. Also the non-diffractive sample is corrected for the con- 
tamination from diffractive events with high dissociated mass; 

• limited track finding efficiency and acceptance of the track selection cuts, as well as the 
limited resolution in momentum and angle; 

• loss of tracks due to secondary interactions and contamination from secondary tracks; 

• decays of charged pions and kaons, photon conversions and decays of lambdas and neutral 
kaons. Thus, in the final spectra the charged kaons appear, while the decay products of 
neutral kaons and lambdas do not. For all the other strange and charmed states, the 
decay products were included. 

The validity of our acceptance correction method relies on the correct simulation of the de- 
scribed effects in the Monte Carlo program. The possible discrepancies between reality and 
Monte Carlo simulation were analyzed and the estimation of the effect on the final distribu- 
tions was included in the systematic uncertainty, as described in the following section. 

8 Systematic effects 

One of the potential sources of systematic inaccuracy is the tracking system and its simulation 
in the Monte Carlo events used for the acceptance correction. Using an alternative simulation 
code with artificially degraded tracking performance we verified that the efficiency to find a 
track which fulfills all the selection cuts is known with an accuracy of about 10%. The error 
due to an imprecise description of the momentum resolution at high px is negligible compared 
to the statistical precision of the data. We also verified that the final spectra would not change 
significantly if the tracking resolution at high px had non-gaussian tails at the level of 10% or 
if the measured momentum was systematically shifted from the true value by the momentum 
resolution. 

Another source of systematic uncertainty is the Monte Carlo simulation of the trigger response. 
We verified that even a very large (20%) inaccuracy of the BCAL energy threshold would not 
produce a statistically significant effect. An incorrect RCAL trigger simulation would change 
the number of events observed, but would not affect the final pt spectra since it is normalized 
to the number of events. The correlation between the RCAL energy and the pt of tracks is 
very small. 

To evaluate the model dependence we repeated the calculation of the correction factors using 
an alternative set of Monte Carlo programs (see section ^) and compared the results with the 
original ones. The differences between the obtained factors varied between 5% for the high 
mass diffractive sample and 11% for the non-diffractive one. The sensitivity of the result to the 
assumed relative cross sections of the physics processes was checked by varying the subprocess 
ratios within the error limits given in |J. The effect was at most 3%. 

All the above effects were combined in quadrature, resulting in an overall systematic uncertainty 
of the charged particle rates as follows: 15% in the non-diffractive sample, 15% in the (Mx) = 
5GeV diffractive sample and 9% in the (Mx) = lOGeV diffractive sample. All these systematic 
errors are independent of px- 
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9 Results 



The double differential rate of charged particle production in an event of a given type is calcu- 
lated as the number of charged particles AiV produced within Arj and Ap T in N ev events as a 
function of Pt- 

1 d 2 N 1 1 AiV 
N ev dp T di] N ev 2p T Ar] Ap T ' 

The charged particle transverse momentum spectrum was derived from the transverse momen- 
tum distribution of observed tracks normalized to the number of data events by means of the 
correction factor described in section |7|. The resulting charged particle production rates in 
diffractive and non-diffractive events are presented in Fig. [I] and listed in Tables [I], |^ and |^. 
In the figure the inner error bars indicate the statistical error. Quadratically combined sta- 
tistical and systematic uncertainties are shown as the outer error bars. The (Mx) = 5 GeV 
diffractive spectrum extends to px = 1-75 GeV and the {Mx) = 10 GeV distribution extends to 
Pt = 2.5 GeV. The non-diffractive distribution falls steeply in the low pt region but lies above 
the exponential fit at higher p? values. The measurements extend to pt = 8 GeV. 
The soft interactions of hadrons can be successfully described by thermodynamic models that 
predict a steep fall of the transverse momentum spectra that can be approximated with the 



exponential form | 21fl : 



d 2 N 



N ev d PT d v - eXp{a ' b -^ + ml) (1) 

where m n is the pion mass. The results of the fits of this function to ZEUS data in the interval 
0.3 < pt < 1.2 GeV are also shown as the full line in Fig. [I]. The resulting values of the 
exponential slope b are listed in Table |J The systematic errors were estimated by varying the 
relative inclusive cross sections within the systematic error limits (see section |8|) and by varying 
the upper boundary of the fitted interval from p T = 1.0 GeV to 1.4 GeV. 

In Fig. ^| we present a comparison of the b parameter resulting from the fits of (|]) to proton- 
proton and proton-antiproton data as a function of the cm. energy. The slope of the ZEUS 
non-diffractive spectrum agrees with the data from hadron-hadron scattering at an energy close 
to the ZEUS photon-proton cm. energy. The diffractive slopes agree better with the hadronic 
data corresponding to a lower energy. In Fig. ^ the ZEUS diffractive points are plotted at 5 GeV 
and 10 GeV, the values of the invariant mass of the dissociated photon. A similar behaviour 
has been observed for the diffractive dissociation of protons, i.e. the scale of the fragmentation 



of the excited system is related to the invariant mass rather than to the total cm. energy |22 
The dashed line in Fig. ^ is a parabola in log(s) and was fitted to all the hadron-hadron points 
to indicate the trend of the data. As one can see, our photoproduction results are consistent 
with the hadronic data. 

The non-diffractive spectrum in Fig. |l] clearly departs from the exponential shape at high px 
values. Such a behaviour is expected from the contribution of the hard scattering of partonic 
constituents of the colliding particles, a process that can be described in the framework of 
perturbative QCD. It results in a high pt behaviour of the inclusive spectrum that can be 
approximated by a power law formula: 

1 d2N A. ( i+^r (2) 



N ev dp T dr] pto' 

where A, pto and n are parameters determined from the data. The fit to the ZEUS points in 
the region of pt > 12 GeV gives a good description of the data and results in the parameter 
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values pto = 0.54 GeV, n = 7.25 and A = 394 GeV -2 . The statistical precision of these 
numbers is described by the covariance matrix shown in Table [|. The fitted function is shown 
in Fig. [I] as the dotted line. 

In Fig. pi the ZEUS data are presented together with the results of a similar measurement from 



the HI collaboration at (W) = 200 GeV [|29| and the data from the WA69 photoproduction 
experiment at a cm. energy of (W) = 18 GeV ||. For the purpose of this comparison, the 
inclusive cross sections published by those experiments were divided by the corresponding total 
photoproduction cross sections |3(| |31|. Our results are in agreement with the HI data. The 
comparison with the WA69 data shows that the transverse momentum spectrum becomes harder 
as the energy of the 7J9 collision increases. Figure |3| also shows the functional fits of the form 



(g) to pp data from UAl and CDF at various cm. energies |27|, [28[. Since the fits correspond 
to inclusive cross sections published by these experiments, they have been divided by the cross 
section values used by these experiments for the absolute normalization of their data. The 
inclusive pr distribution from our photoproduction data is clearly harder than the distribution 
for pp interactions at a similar cm. energy and in fact is similar to pp at = 900 GeV. 
This comparison indicates that in spite of the apparent similarity in the low px region between 
photoproduction and proton-antiproton collisions at a similar cm. energy, the two reactions are 
different in the hard regime. There are many possible reasons for this behaviour. Firstly, both 
of the pp experiments used for the comparison measured the central rapidity region (|?7| < 2.5 
for UAl and \r]\ < 1 for CDF), while our data correspond to 0.8 < r\ cm . < 3.4. Secondly, 
according to VDM, the bulk of the ^p collisions can be approximated as an interaction of a 
vector meson V with the proton. The pt spectrum of Vp collisions may be harder than pp at a 
similar cm. energy, since the parton momenta of quarks in mesons are on average larger than in 
baryons. Thirdly, in the picture where the photon consists of a resolved part and a direct part, 
both the anomalous component of the resolved photon and the direct photon become significant 
at high pt and make the observed spectrum harder compared to that of Vp reactions. 
Figure |] shows the comparison of our non-diffr active data with the theoretical prediction 
obtained recently from NLO QCD calculations [32]. The charged particle production rates in a 
non-diffractive event were converted to inclusive non- diffract ive cross sections by multiplying 
by the non-diffractive photoproduction cross section of a n d{.lP — ► X) — 91 ± 11 fih 0. The 
theoretical calculations relied on the GRV parametrisation of the parton densities in the photon 
and on the CTEQ2M parametrisation for partons in the proton p3|. The NLO fragmentation 
functions describing the relation between the hadronic final state and the partonic level were 
derived from the e + e~ data [|33j. The calculation depends strongly on the parton densities in 
the proton and in the photon, yielding a spread in the predictions of up to 30% due to the 
former and 20% due to the latter. The factorization scales of the incoming and outgoing parton 
lines, as well as the renormalization scale, were set to pt- The uncertainty due to the ambiguity 
of this choice was estimated by changing all three scales up and down by a factor of 2. The 
estimates of the theoretical errors were added in quadrature and indicated in Fig. |asa shaded 
band. The theoretical calculation is in good agreement with the ZEUS data. 



10 Conclusions 

We have measured the inclusive transverse momentum spectra of charged particles in diffractive 
and non-diffractive photoproduction events with the ZEUS detector. The inclusive transverse 
momentum spectra fall exponentially in the low px region, with a slope that increases slightly 
going from the non-diffractive to the diffractive collisions with the lowest Mx- The diffractive 



8 



slopes are consistent with hadronic data at a cm. energy equal to the invariant mass of the 
diffractive system. The non-diffractive low pt slope is consistent with the result from pp 
at a similar cm. energy but displays a high p? tail clearly departing from the exponential 
shape. Compared to photoproduction data at a lower cm. energy we observe a hardening 
of the transverse momentum spectrum as the collision energy increases. The shape of our px 
distribution is comparable to that of pp interactions at y/s = 900 GeV. The results from a 
NLO QCD calculation agree with the measured cross sections for inclusive charged particle 
production. 
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MGeV] 



cFN 



N ev dp^dr) 



GeV~ 



0.30- 
0.40- 
0.50- 
0.60- 
0.70- 
0.80- 
0.90- 
1.00- 
1.10- 
1.20- 
1.30- 
1.40- 
1.50- 
1.60- 
1.70- 
1.80- 
1.90- 
2.00- 
2.14- 
2.29- 
2.43- 
2.57- 
2.71- 
2.86- 
3.00- 
3.25- 
3.50- 
3.75- 
4.00- 
4.50- 
5.00- 
5.50- 
6.00- 
7.00- 



0.40 
0.50 
0.60 
0.70 
0.80 
0.90 
1.00 
1.10 
1.20 
1.30 
1.40 
1.50 
1.60 
1.70 
1.80 
1.90 
2.00 
2.14 
2.29 
2.43 
2.57 
2.71 
2.86 
3.00 
3.25 
3.50 
3.75 
4.00 
4.50 
5.00 
5.50 
6.00 
7.00 
8.00 



4.98 
2.99 
1.78 
1.09 
0.641 
0.420 
0.259 
0.164 
0.107 
0.0764 
0.0513 
0.0329 
0.0242 
0.0175 
0.0133 
0.0082 
0.00615 
0.00454 
0.00360 
0.00215 
0.00166 
0.00126 
0.00098 
0.000625 
0.000456 
0.000252 
0.000147 
0.000094 
0.000067 
0.0000301 
0.0000151 
0.0000082 
0.0000038 
0.0000014 



VstaA GeV" 



&syst[ GeV" 



0.05 


0.74 


0.03 


0.44 


0.02 


0.26 


0.01 


0.16 


0.012 


0.096 


0.010 


0.063 


0.007 


0.038 


0.005 


0.024 


0.004 


0.016 


0.0034 


0.0114 


0.0017 


0.0077 


0.0012 


0.0049 


0.0010 


0.0036 


0.0008 


0.0026 


0.0006 


0.0020 


0.0005 


0.0012 


0.00038 


0.00092 


0.00028 


0.00068 


0.00024 


0.00054 


0.00017 


0.00032 


0.00013 


0.00025 


0.00012 


0.00018 


0.00010 


0.00015 


000071 

V 7 • V 7 V 7 V 7 V 7 1 -L 


000093 

V 7 • V 7 V 7 V 7 V 7 7 ' J 


0.000048 


0.000068 


0.000031 


0.000037 


0.000020 


0.000022 


0.000012 


0.000014 


0.000008 


0.000010 


0.0000045 


0.0000045 


0.0000029 


0.0000023 


0.0000021 


0.0000012 


0.0000009 


0.0000006 


0.0000005 


0.0000002 



Table 1: The rate of charged particle production in an average non-diffractive event. The data 
correspond to —1.2 < r] < 1.4. The a stat and a syst denote the statistical and systematic errors. 
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i .^V[GeV- 2 " 




<j,tnt\ GeV~ 2 l 


a syst [GeV T 




0.30- 0.40 


1.63 


0.06 


0.24 


0.40- 0.50 


1.02 


0.04 


0.15 


0.50- 0.60 


0.559 


0.028 


0.083 


0.60- 0.70 


0.308 


0.019 


0.046 


0.70- 0.80 


0.165 


0.013 


0.024 


n fin n on 
U.oU— u.yu 


0.088 


U.U1 1 


0.013 


0.90- 1.00 


0.0479 


0.0059 


0.0071 


1.00- 1.10 


0.0312 


0.0052 


0.0046 


1.10- 1.20 


0.0196 


0.0042 


0.0029 


1.20- 1.35 


0.0100 


0.0018 


0.0015 


1.35- 1.50 


0.00304 


0.00087 


0.00045 


1.50- 1.75 


0.00153 


0.00052 


0.00023 



Table 2: The rate of charged particle production in an average event with a diffractively disso- 
ciated photon state of a mass (M x ) = 5 GeV. The data correspond to —1.2 < r\ < 1.4. The 
&stat an d o S yst denote the statistical and systematic errors. 



MGeV] 



N e . 



dp^drj 



GeV 



-2 



o- stat \GeV 



-2 



a syst [ GeV 



-21 



0.30- 
0.40- 
0.50- 
0.60- 
0.70- 
0.80- 
0.90- 
1.00- 
1.10- 
1.20- 
1.35- 
1.50- 
1.75- 
2.00- 



0.40 
0.50 
0.60 
0.70 
0.80 
0.90 
1.00 
1.10 
1.20 
1.35 
1.50 
1.75 
2.00 
2.50 



3.87 
2.32 
1.46 
0.803 
0.485 
0.288 
0.176 
0.109 
0.0732 
0.0294 
0.0186 
0.0086 
0.00260 
0.00076 



0.10 


0.34 


0.06 


0.20 


0.04 


0.13 


0.033 


0.072 


0.023 


0.043 


0.017 


0.025 


0.012 


0.015 


0.009 


0.009 


0.0075 


0.0065 


0.0035 


0.0026 


0.0028 


0.0016 


0.0014 


0.0008 


0.00066 


0.00023 


0.00023 


0.00007 



Table 3: The rate of charged particle production in an average event with a diffractively disso- 
ciated photon state of a mass (M x ) = 10 GeV. The data correspond to —1.2 < r\ < 1.4. The 
&stat an d a syst denote the statistical and systematic errors. 
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sample 


b[ GeV" 1 ] 


&stat(b) 


&syst(b) 


a 


v s tat{a) 


cov(a, b) 


non-diffr active 


4.94 


0.09 


0.19 


3.39 


0.09 


-0.011 


diff (M x ) = 5 GeV 


5.91 


0.17 


0.19 


2.78 


0.10 


-0.016 


diff (M x ) = 10 GeV 


5.28 


0.10 


0.17 


3.34 


0.06 


-0.006 



Table 4: The values of the parameters resulting from the fits of equation to ZEUS data in 
the interval 0.3 < px < 1-2 GeV. The a s t a t and o- syst indicate the statistical and systematic 
errors. 





Pto 


n 


A 


Pto 


0.32 • 10- 3 


0.48 ■ 10- 3 


-0.10- 10 1 


n 




0.12 ■ 10- 2 


-0.12 ■ 10 1 


A 






0.32 ■ 10 4 



Table 5: The covariance matrix corresponding to the fit of equation (^) to the non-diffr active 
data for p T > 1.2 GeV. 
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Figure 1: Inclusive transverse momentum distributions of charged particles in photoproduction 
events at (W) = 180 GeV averaged over the pseudorapidity interval of —1.2 < rj < 1.4. The 
inner error bars indicate the statistical errors and the outer ones represent the quadratic sum 
of the statistical and systematic errors. Solid lines indicate fits of equation (j7]j to the data 
in the region of px < 1-2 GeV. The dotted line shows a power law formulae (TJj fitted to the 
non- diff r active data forpx > 1-2 GeV. For the sake of clarity the diffractive points are shifted 
down by two orders of magnitude. 
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Figure 2: Inverse slope of the exponential fit of the form (^) vs. the cm. energy for different 
experiments. The fixed target data were taken from Jjffij J^/ /p3|/, ISR - fftBJ, UA1- The 
ZEUS non-diffractive point is placed at the photon-proton cm. energy. The diffractive points 
are plotted at the energies corresponding to the mean value of the invariant mass (Mx) ■ The 
error bars indicate the quadratic sum of the statistical and systematic errors. The dashed line 
is a parabola in log(s) and was fitted to all the hadron-hadron points to indicate the trend of 
the data. 
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Figure 3: Comparison of ZEUS non-diffractive transverse momentum spectrum with the data 
from HI OMEGA /§/, UA1 Jjgjj and CDF The inner error bars indicate the statistical 
errors and the outer ones represent the quadratic sum of the statistical and systematic errors. 
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Figure 4: Comparison of the ZEUS inclusive cross sections for non-diffractive photoproduction 
at (W 1P ) = 180 GeV with the NLO QCD calculation results from fffi^ . The inner error bars on 
the data points indicate the statistical errors and the outer ones represent the quadratic sum 
of the statistical and systematic errors. The shaded band corresponds to the uncertainty of the 
theoretical calculation. 
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